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Abstract 

The second-order QCD matrix elements give a very good agreement with 
experimental data on the angular distributions of the four-jet events in e^e“ 
collisions at the resonance energy. Unfortunately the description of the 
sub-jet structure is quite poor. The alternative approach, parton showers, 
gives a good description of the sub-jet structure but is worse than matrix 
elements when it comes to the angular distributions. Here is presented a 
hybrid between the matrix elements and the parton showers that is intended 
to combine the best of the two approaches. 


1 Introduction 


According to the viewpoint of modern particle physics all matter in the universe is com¬ 
posed of the particles listed in the table below. 


Charge 

Quarks 

2/3 

-1/3 

Leptons 

-1 

0 

Generation 

hrst 

u (5 ■ 10-3) 

d (7-10-3) 

e (0.5 ■ 10-3) 

( 0 ) 

second 

c (1.5) 

s ( 0 . 2 ) 

p (0.106) 

( 0 ) 

third 

t (170) 

b(5) 

r (1.784) 

Ur ( 0 ) 


The approximate particle masses in the parentheses are given in GeV. All these parti¬ 
cles have spin 1/2 and are therefore fermions. This means that the total wave-function 
'l/(gi,..., gAr) of a iV particle system is antisymmetric under the permutation of two 
identical particles i.e. 'h(..., ..., • • •) = —d'(- ■ ■ ,(lj ■ ■ ■ ,(li, ■ ■ ■) ■ Thus two identical 

particles i and j can’t have identical quantum numbers g* = qj because this would imply 
that T = 0. This is the so-called Pauli principle that determines the structure of the pe¬ 
riodic system. The particles in the second and third generations have the same properties 
as the corresponding particles in the hrst generation except for the masses that increases 
with the generation number. It is only the particles in the hrst generation that is needed 
to build the matter on earth. 

There are four known fundamental interactions between the matter particles. These 
are the electromagnetic, the weak, the strong and the gravitational interactions. For 
the hrst three interactions we have quantum mechanical theories possessing local gauge 
symmetries, in which the forces are mediated by exchange of virtual particles. The fourth 
interaction, gravitation, is currently described by Einstein’s general theory of relativity, 
which is a classical theory. One of the main objectives of current theoretical research is 
to hnd the correct quantum mechanical description of the gravitation. Fortunately the 
gravitational interactions between the particles studied in particle physics experiments 
are so weak that they can be leaved out of account. The interactions, their mediating 
particles and the matter they affect are listed in the table below. 


Interaction 

Particle 

Affects 

electromagnetic 

7 

all except z/g , 

weak 

w+ w- 

all particles 

strong 

gi * = 1 • 

.. 8 all quarks 

gravitation 

graviton 

all particles 


All these particles have integer spin and are therefore bosons. The 7 particle which medi¬ 
ates the electro-magnetic interaction is a spin 1 particle. The virtual 7 particle is massive 
and can therefore exist in the tree different spin-states — 1,0 and -|- 1 , this corresponds to 
the fact that the classical E and B helds can be described by three independent parame¬ 
ters. Similar arguments lead to the conclusion that the graviton, if it exists, should be a 
spin 2 particle. The electromagnetic interaction is described by QED (Quantum Electro 
Dynamics) a theory whose predictions agrees extremely well with experimental results. 
This is in part due to the smallness of a = e^/dvreohc in which the perturbative solutions 
are expanded. The weak interaction is described by a theory that combines the electro¬ 
magnetic and the weak interaction into a theory of the so called electro-weak interaction. 
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In the electro-weak theory the fundamental particles that mediate the interactions are 
the B, W°, W"*" and W“ particles. They are all massless, a condition that is necessary for 
a consistent theory. The W’*' and W“ particles acquire their masses through interaction 
with the Higgs field that hlls all of space. This is analogous to the repeated absorption 
and emission of photons e.g. in glass, which slows down the overall velocity of the pho¬ 
tons and thus gives them an effective mass that is larger than zero. The same Higgs held 
polarizes the superposition of B and W° into the physical states 7 and TP, where the 7 
remains massless and the Z° acquires a mass. 

The theory of the strong interaction is described by QCD (Quantum Chromo Dynam¬ 
ics) a gauge theory based on the symmetry group SU(3). The group SU(3) could in matrix 
language be described as the set of unitary 3x3 matrices whose determinants equals unity. 
Such a matrix has 8 independent parameters that corresponds to the eight gluons that 
mediates the interaction. The quarks, which are the only particles to be affected by the 
strong interaction, carry a strong charge. This comes in three different types, whimsically 
labeled R(ed) G(reen) and B(lue). The antiquarks carries the corresponding anticharges 
labeled R, G and B. 

Quarks never exists freely, but are always combined into color singlet states called 
hadrons. There are three fundamental ways to build a hadron: 

1. equal combinations of colors RGB 

2. equal combinations of anticolors RGB 

3. equal combinations of a color and its anticolor RR -|- GG -|- BB 

No 1 corresponds to baryons, no 2 to antibaryons and no 3 to mesons and antimesons. 
The quarks and gluons are collectively called partons, this term stems from the days 
before the quark model was well established. QGD is a non-Abelian held theory; this 
corresponds to the fact that the gluons carry color charge themselves. The gluons can 
thus interact with each other. This self-interaction give rise to the conhnement of quarks 
into hadrons. The held lines of the color held between two quarks forms a narrow tube 
when the quarks separation increases, instead of spreading out as the electromagnetic 
held lines. The energy in the color held is thus proportional to the distance between the 
quarks. When the potential energy is big enough a quark antiquark pair is created. The 
consequence is that there can never exist a free color charged quark. Another consequence 
of the self-interactions is that the equations of QGD can only be solved perturbatively for 
short distances, i.e. high energies. For larger distances the perturbative approach fails. 
Since the equations of QGD are too complex to be solved directly one must use some kind 
of simplihed picture. One such approximation is the Lund model. 

This article is focused on the QGD treatment of the e’''e“ —hadrons event. The article 
is organized as follows: Section 2 presents the two standard methods used to describe the 
perturbative phase of QGD, i.e. matrix elements and parton showers. It also contains a 
description of fragmentation of quarks into hadrons and the jets that these hadrons form. 
Section 3 describes the hybrid of matrix elements and parton showers, the main result 
of my diploma work. Section 4 shows the theoretical predictions of the hybrid compared 
with those of the matrix elements and the parton showers. 
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Figure 1: e’''e 


qq 
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Existing descriptions of the e+e 


hadrons event 


One experimental way to study QCD is the collision of e"*" with e“. This method has the 
advantage that the process is very clean in the sense that both e+ and e“ are considered 
to be fundamental point-like particles. The process could to hrst order in as be described 
by the Feynman diagram in hgure The zeroth order cross section cxo is approximately 
given by 


(To ~ 


/A B 

_ m|o)2 m|oF|o 


-|- interference 


( 3 ) 


where p is the four momentum of either 7 or The zeroth order cross section contains 
no QCD corrections except for a factor 3 counting the number of possible quark colors. 
The hrst term corresponds to the 7 channel while the second term corresponds to the Z° 
channel. There is also a 7 /Z° interference term between the two channels. The factors A 
and B depends on the electro-magnetic and the weak coupling constants. The quantity 
p^ equals the total energy in the center of momentum frame that coincides with the rest 
frame of the detector. The cross section has two peaks at = 0 GeV and m^o = 91.2 
GeV. Experiments designed for detailed studies of QGD are therefore preferably run at 
Ecm = 91.2 GeV in order to obtain maximum statistics. 


2.1 Matrix elements 

The perturbative solution to the QGD equations is expanded in terms of Og, the strong 
coupling constant. The matrix elements are the terms of this expansion. They can 
be represented by Feynman diagrams. Each three-parton vertex involving the strong 
interaction is assigned a factor The four-gluon vertex has a factor as but it does not 
appear until in higher orders than considered here. Because the probability amplitudes 
are squared to obtain the cross section the number of strong interaction vertices equals 
the order in (Ug. The strong coupling constant is, despite its name, a function of 
where cm refers to the center of momentum frame of the whole event. To hrst order ctg is 
given by 0 

^ (33 - 2uf) log(Q7A2) • 

The integer Uf is the number of quark havors available at the given energy, typically 4 or 
5. As seen from expression (^) as{Q‘^) decreases with increasing energies and therefore 
becomes small for short-distance interactions. We say that QGD is asymptotically free. 
The parameter A is a free parameter that has to be determined by experiments. According 
to experiments A lies in the interval 0.2 to 0.3 GeV . 
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The matrix elements are in this article only showed explicitly to hrst order in Og 
becanse the second order matrix elements 0 are far too complex. The amplitnde for a 
3-parton event is to first order given by the Feynman diagrams in hgnre |^. The cross 


e 



q 

g 

q 


e 



q 

g 

q 


Fignre 2: e'''e —> qgq 


section ai is given explicitly as 

1 dcTi _ as 4 xl + xl _ 2Ei . x 

(To dxi dx2 27r 3 (1 — Xi)(l — X 2 ) * T^cm ’ ’ 

where 1,2,3 corresponds to q,q,g and the energies Ei are given in the center of momentum 
frame. The masses of the partons are assumed to be 0. It is then easy to show that 
0 < Xj < 1. Expression (|^) is divergent in the region where at least one of xi or X 2 
approaches 1. The case where only one of them approaches 1 is called the collinear 
singularity, because the other two parton momenta become collinear. The other case 
when both approaches 1 is called the soft singularity because the gluon energy becomes 
0. The total first order cross section atot is given by 

(Ttot “ <^0 = <^0 + C’’! + <^virt (6) 

where is the first order virtual corrections created by the interference between the 
diagrams in figure and The above mentioned divergences in ai are cancelled by 



Figure 3: first order virtual corrections 


divergences in a^irt thus rendering atot finite. The divergences in at is therefore no real 
analytical problem but the use of Monte Carlo methods requires 0 < at/atot <1- To 
handle this we introduce cuts in the phase space by demanding that yjk > 2/min • The 
variable yjk is given by y^k = {pj + Pj )"^/= 1 — x* where indices i,j,k are all different. 

By plotting the expression 




' Vij >ymin 


1 d(Ti 
(To dxi dx2 


dxi dx2 


^ 1 

27r 3 


Vij >2/min 2/23 2/13 


(7) 
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as a function of j/min with ag ~ 0.15 one notices that must be at least 0.01 to ensure 
that (Ji/(Ttot < 1 • This corresponds to a minimum mass 


mmin 




minT/cm 


( 8 ) 


At the resonance we obtain mmin — 9 GeV. The running of Og sets a limit for the 
perturbative approach at ^ 1 GeV. The result is that the gluons in the region 1 to 9 
GeV are missed. As seen from expression (|]) the region of discarded gluons increases as 
a linear function of the energy T^cm- This is a severe shortcoming of the matrix element 
approach. 

The amplitudes for a 4-parton event to second order in ttg is given by the Feynman 
diagrams listed in £g The diagrams can be divided into three categories: double- 

1 2 3 





Figure 4: 4-parton histories to second order in Og 


bremsstrahlung (1-3), triple-gluon vertex (4-5) and secondary qq-production (6-7). The 
cross sections for reaction e’''e“ —*• 7 */Z° —qqgg and reaction e+e“ —>• 7 */Z° —> qqqq are 
given by 0 


‘^qqgg(2/u) 

<7qqqq(l/p) 


C^O 


C^O 


Cf^Aivij) + (g| - IcfCa) B(yij)+C^C^C(yij) 
CrTrUiDlyij) + (Cp — -CpCa) E{yij) 


(9) 

( 10 ) 


where Uf is the number of active quark flavors and i/ij = is the normalized two- 

body invariant mass with indices i and j running over the four partons. The functions 
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Aijjij)... Eijjij) are group independent and contains the full kinematical dependence of 
the cross section. A set of hve independent yij are needed to specify the kinematical con- 
hguration of a 4-parton event. The so-called Casimir factors C-p, Ca, Tp are respectively 
a measure of the coupling strengths of the reactions q —>• qg, g —> gg and g —> qq. Thus 
the hrst term in expression correspond to diagrams 1 to 3 in hgure ^ and the third 
term correspond to diagrams 4 and 5. The second term in expressions @ and O are 
formed by the quantum mechanical interference. The Casimir factors Cp, Ca, Tp are 
given directly by the symmetry group. For QCD the group is SU(3) and the Casimir 
factors have the values Cp = 4/3, Ca = 3 and Tp = 1/2. 

2.2 Parton showers 

A parton shower is a kind of semi-classical approximation of parton events in the sense 
that every parton in the event has a given history with a specihed four-momentum. Under 
the simulation each particle is allowed to branch into two new particles. The three basic 
branchings are shown in hgure By iterating these basic branches a hnal state with an 

q ^ qg g ^ gg g ^ qq 




Figure 5: The basic branchings in a parton shower 


arbitrary number of partons may be constructed. The differential probabilities for each 
branching are given by the so-called Altarelli-Parisi equations of evolution |^: 


dFa^bc — 


27t 




( 11 ) 


The splitting kernels Pa^bc{z) are given by the following expressions depending on the 
type of branching: 


-^q^qg(^) 

^g-gg(^) 

^g-qq(^) 


4 

3 

3 


l + z"^ 
1-z ’ 


5 


Hf 


(^2 + (1 - 


( 12 ) 


The splitting kernels can be obtained from the 3-parton matrix elements. By dehning 
z = Xi expression turns into: 


dR 


da 


Ofo 


^2 I ^2 
Jy 2 *^2 


q^qg 


(To 27r 3 (1 — xi)(l 
^ 4 z^ + (1- yisf 
27r 3 (1 - z)yi3 


X 2 ) 
dz dyi 3 


dxi dx2 
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a. 


4 dmq 


+ 1 


27r 3 1 


dz 


a 


„ dm? 
s q p 

f 


27r 


q^qg 


»dj. 


( 13 ) 


In the last step comparison with expression (p!TD gives the expression for Pq^qg( 2 :). The 
expressions Pg^gg and Pg^qq could be obtained in a similar way. As seen from the approx¬ 
imation in expression ( 0 ) the parton shower only gives a good description in the limit 
1/13 ^ 0, i.e. when the partons produced in the branching becomes collinear. Probabilities 
larger than unity in the collinear region correspond to possibility of a quark emitting more 
than one gluon. Parton showers are therefore complementary to the matrix elements in 
the sense that they are well behaved below the i/min cut at 9 GeV. The exact definition 
of z and varies amongst different algorithms but generally z is some variant of energy 
sharing between the two daughters and is a function of and z. Furthermore the 
expression dmq/mq in equation (pd^ can be replaced by [z))/{ml^f {z)) where f{z) 

is an arbitrary well-behaved function of z. 

The probability that a parton does not branch between an initial mass m and a 
minimum value mmin is given by the integration and exponentiation of expression ([TT|) : 


Saijn) = exp 


dm'2 as(Q‘^) 

Zmin(m') 27r 


/m^ . m 

min 


dzPa^hc{z) 


(14) 


where is a small cutoff mass that is used to regularize the collinear divergences in 
z and the infrared divergences in m^. Expression ([T^) is called the Sudakov form factor. 
An analogue to expression (|T^ is the formula 


which gives the fraction of radioactive nucleus still remaining after a time t. The proba¬ 
bility that a parton with a maximum mass rumax will branch in the interval [m^, 
is given by 

Pa(mLx,"i') = Sa{mU-^ 

A branching with a specific m is then selected using a Monte Carlo method. A random 
number R is chosen uniformly in the interval (0,1) and the equation 


Saim^) = (16) 

is solved for m. Unfortunately is it often impossible to solve equation (0) analytically 
but appropriate numerical methods exists. Whenever the lower limit m = is reached 
the particle is put on its mass-shell and the evolution of the particle is terminated. 

An independent evolution of each parton in the shower overestimates the total amount 
of evolution. The coherence effect can to some extent be included by imposing angular 
ordering, i.e. by requiring that the angle between two daughter partons decreases as one 
goes from the initial to the final partons. A consequence of the angular ordering is a slower 
multiplicity growth of partons and a depletion of parton production at small energies. 
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As previously mentioned the definition of and 2 ; in the parton shower varies between 
the different parton shower algorithms. In JETSET |^, the event generator that I’ve 
been using, the argument in as is given by 

= z{l - z)ml (17) 

where Pt is the component of the daughters momentum transverse to the momentum of 
the parent. The definition of is motivated by the attempt to include coherence effects. 
To choose a definition for z we first look at the system in the rest-frame of parton a and 
with the z-axis collinear to the velocity of a in the parton shower CM-frame, see figure |^. 
If we now translate to the parton showers CM-frame it is straightforward to show that 



Z-axis 


Figure 6: The branching a —> 6c in the rest frame of a. 


Eh 1 + pcos9 (Pa)^ 

~ Ea~ 2 ^ 


(18) 


under the assumption that b and c are massless. In the infinite-momentum frame, (3^1, 
this reduces to 2 ; = (1 -|- cos6')/2. If we now decide that relation (11^) should hold even if 
b and c are massive then, after some calculations, z turns out to be given by: 


_ ml Eh ml — X +ml — ml 

^ ~ 2A 

with A = {ml — ml — mlY — Ami • 


(19) 


2.3 Fragmentation 

At large times, when the partons produced by a matrix element or a parton shower 
become more separated, the strong coupling constant Og increases to the limit where 
perturbative calculations become impossible. Because of the increasing strength of the 
strong interaction all color-charged partons are forced to combine into colorless states, the 
so-called hadrons. It is this process of free quarks and gluons turning into hadrons that is 
called fragmentation. The process is not yet understood in the context of the fundamental 
QCD equations, so some kind of phenomenological method is needed. These methods can 
be divided into 3 categories: independent fragmentation, string fragmentation and cluster 
fragmentation. One example of string fragmentation is the ‘Lund model’. To explain the 
string model we first look at how it works in the simplest case with a qq color-singlet. 
Calculations with lattice QCD indicate that the energy in the color-field is proportional 








to the distance between the partons. This leads to the picture of a color flux tube that 
is stretched between the two particles. The flux tube is uniform along its length and has 
a transverse size of roughly 1 fm. When the distance between the endpoints q and q 
increases the energy in the tube becomes large enough to create two new partons q'q'. 
The result is two new color-singlets qq' and q'q each with its own flux tube. The process 
is then repeated until only on-mass-shell hadrons remains. 

In the more complex situation with the initial conhguration qgq, shown in hgure |^, 
one draws a tube from the quark to the gluon and then to the anti-quark. The gluon 


g 



Figure 7: String drawing for a qgq event 

must be placed inside the string because the gluon has a color-charge given by a color 
and a anti-color. The charges of the string q - g - q is of the form R - RG - G. The 
sub-strings q - g and g - q are then evolving as in the qq example except for one meson 
that is produced at the gluon-corner consisting of one parton from each sub-string. When 
the partonic state contains 2n quarks then n strings are formed, each containing a quark, 
an antiquark and an arbitrary number of gluons. 

The procedure above only produces mesons. One easy way to obtain the production 
of baryons is to treat a diquark in a color antitriplet state as an ordinary antiquark. The 
strings can then either break by quark-antiquark or antidiquark-diquark pair production. 
An alternative model is the ‘popcorn’ model in which the antidiquark-diquark pair pro¬ 
duction is replaced by two consecutive quark-antiquark pair productions. The ‘popcorn’ 
model gives a less strong correlation between the color and momentum space of the baryon 
and the antibaryon coming from the same pair production. 

Many of the particles produced in the fragmentation process are unstable particles 
with very short lifetimes. The simulation of the decay process is governed by tables 
over the possible decay modes and their corresponding branching ratios. Normally it is 
assumed that decay products are distributed according to phase space, i.e. there is no 
dynamics involved in their relative distribution. 

2.4 Jets and their angular distributions 

The angular distribution of the particles detected in an e’''e“ collision is not uniform. 
Instead the particles are clustered together in so-called jets. The directions of the jets 
are strongly correlated to the directions of the hardest partons in the perturbative phase. 
Since there is a continuum of differently separated parton emissions there is no unique 
jet dehnition. In order to determine the number of jets and their corresponding four- 
momenta some kind of cluster algorithm has to be used on the hnal particles. One type of 
cluster algorithm are those based on binary joining. In these algorithms all particles are 
considered to be separate clusters at the start. Then a distance measure is calculated 
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for every pair of clusters. If the smallest dij is smaller than some given djom then the 
corresponding clusters are joined into a new cluster. The procedure is then repeated for 
the remaining clusters until all clusters are separated by a distance greater than djom- 


The final clusters are then the seeked-for jets. One possible definition of dij used in the 
JADE scheme IHl is 


dij — 


2EiEj{l — cos 9 ij) 


El 


( 20 ) 


Here E^^s is the total visible energy of the event. In an event simulation and in an ideal 
experiment E^ig equals Ecm- The dimensionless nature of dij in expression (^0|) makes 
it suitable for comparing results at different c.m. energies. When given the jets from a 
4-parton event there is no unique pairing of partons with jets. The quark-jets are harder 
than the gluon-jets on the average, but the difference smears out for the individual jets 
under the fragmentation and the clustering process. With the four jets enumerated in 
decreasing energy order the hrst two jets preferentially correspond to the initial qq-par 
while the two last jets preferentially correspond to the emitted gluons or the secondary 
produced qq-par. The uncertainty in the parton-jet pairing prevents one from accessing 
the full hve-dimensional distribution of the 4-parton event. Instead there is need for 
simpler measures that are less sensitive to the above discussed uncertainties. The following 
angular distributions for 4-jet events are examples of such measures: 


1. The Bengtsson-Zerwas correlation |^: 


cosxbz = 


l(Pi X P 2 ) • (P3 X Pi)\ 
\Pi X P2IIP3 X P4I 


( 21 ) 


2. The modified Nachtmann-Reiter variable |§]: 


Icos^nr 


I (pi - P 2 ) ■ (P3 - Pa) I 

\Pl - P2\\P3 - PaI 


( 22 ) 


3. The cosine of the angle between the two jets with lowest energy 0: 


COStt34 — 


P3-P4 

\P3\\P4\ 


(23) 


The three variables are not orthogonal but still sufficiently different to provide comple¬ 
mentary information. The variable cosyez measure the angle between the plane spanned 
by the two most energetic jets and the plane spanned by the two least energetic jets. 
The g —>■ gg vertex in QCD tends to align the two planes, contrary to the g —>• qq ver¬ 
tex that favors a perpendicular relation between the planes. This is due to the different 
spins of the quark and the gluon. The variable could therefore be used to measure the 
difference between QCD and an Abelian variant of QCD that lacks the g —>■ gg vertex. 
The Abelian variant is not a serious competitor to QCD because it doesn’t explain the 
running of as, it is rather a counter-example used to prove the non-Abelian character of 
QCD. The original Nachtmann-Reiter variable measure the angle 0 nr between the two 
vectors pi and p^. The g —>• gg vertex tends to decrease the 6 *nr angle contrary to the 
g ^ qq that tends to increase ^nr. This is again due to the different spins of the quark 
and the gluon. The modihed angle 6 *nr becomes equal to 6 *nr in the limit when parton 
1 and 2 are back-to-back whereas when the angles between parton 1 and 2 or between 


10 







parton 3 and 4 decrease, is roughly equal to xbz- The observable is thus com¬ 
plementary to Xbz- It is used in order to distinguish between two-gluon dual states and 
secondary qq-production. The variable cosa 34 is used in order to distinguish between the 
double-bremsstrahlung and the triple-gluon vertex. 


3 Creating a forced parton shower 


The second-order QCD matrix elements give a very good agreement with experimental 
data on the angular distributions of the four-jet events in e''“e“ collisions at the Z° res¬ 
onance [T^. Unfortunately the description of the subjet structure is quite poor. The 


alternative approach, parton showers, gives a good description of the subjet structure but 
is worse than matrix elements when it comes to the angular distributions. To enhance 
the theoretical accuracy there are different paths to follow. One way is to calculate ma¬ 
trix elements to a higher order. Unfortunately this is impossible in practice because of 
the enormous amount of work needed. Another way is to first use the matrix element 
and then apply a parton shower on the resulting partons. The problem here is that the 
parton shower is normally formulated so that it starts from two initial partons in order 
get all its internal procedures right. For example the angular ordering procedure doesn’t 
have any initial angle to start from. In this report the attempt is to combine the matrix 
element and the parton shower in a more subtle way. The combination is done in the 
sense that a parton-shower is forced to follow a history derived from a matrix-element in 
the first steps of its evolution. The forced parton shower has been implemented as a set 
of Fortran?? subroutines which in turn heavily rely on the event generator JETSET 
?.410. The whole procedure is described in detail in the following text: 


1. The 4-parton matrix element procedure based on the paper p| is executed for 
given values on Bcm and ymm- The result is the four-momentum and the particle 
codes of the final particles in either reaction e’''e“ —> 7 */Z° —>■ qqgg or in reaction 
e’''e“ —*• '-y*/TP qqqq- How these final particles were formed is not known; this 
is due to the underlying quantum mechanical theory that treats the process as a 
superposition of all possible histories. 


2 . 


The Feynman diagrams corresponding to all 4 jet histories of second order in Os 
are shown in fig |^. The enumeration of the diagrams is the same as is used in the 
program source code to label the different histories. The particle codes are examined 
in order to decide if the event, generated in step 1, was a qqgg, event number 1-5, 
or a qqqq event, number 6 -?. When the event type has been narrowed down to one 
of these categories it still remains to select one of the relevant histories. The exact 
choice of history can not be uniquely determined, due to the quantum mechanical 
interference, but a probability can be assigned to each history. The probability for a 
diagram is given by the product of the branching probabilities for the two vertices. 
The branching probabilities are given by expression (pTl) but with the approximation 
(necessary for the matrix-element approach) that as{Q‘^) is assumed to be equal for 
the two vertices. The phase space is chosen as dm'^dz in order to have a consistency 
between the parton-shower algorithm and the matrix-element descriptions, cf. eq. 
(p!3D. For example the relative probability for the history shown in figure is given 
by: 


1 4 1 + 1 g (1 ~ ^56(1 ~ ^se))^ 

mf 31- 2:34 ml ^56(1-^56) 


(24) 
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where 


= pI = {Ps+Pb+Pef , (25) 

ml = pl = {p5+P(>f, 

and z values are evaluated using expression (JT^). Once the probabilities have been 



Figure 8 : 4-parton history no 4 

calculated a specihc history is chosen in one of the following ways: 

• Monte Carlo: a history is chosen at random but according to the relative 
probabilities calculated above. 

• Winner takes all: the history with the largest probability is chosen. 

3. Once a history is chosen the values of z, ml and a </ 3 -angle is recorded for each of 
the two vertices. The two 93 -angles are obtained in the following way. First the 
coordinates are rotated so that the initial quark momentum is in the positive z 
direction, then a second rotation is done so that the the mother of the considered 
vertex has its momentum in the positive z direction. These rotations consist of 
a rotation in the a:|/-plane followed by a rotation in the xz-pla.ne. The need to 
specify a specihc rotation procedure is due to the fact that the orientation of a rigid 
body is determined by three Euler angles. The (/ 9 -angle is then given by the polar 
representation of one of the daughters momentum: 

Px = r sin 6 * cos (/9 (26) 

Py = r sin 6 sin ip 
Pz = r cos 9 

In a q—>qg or a q —0 qg vertex the daughter is chosen to be the quark. 

4. A modihed version of the parton shower procedure in JETSET is executed. Given 
the chosen history the procedure forces the z-value and the mGvalue of the two 
vertices in the corresponding diagram to be as close as possible to the values given 
by the matrix element. The forced ^ value and the z value given by the matrix 
element differs slightly, see hgure This is due to the fact that the hnal partons 
given by the matrix element is on the mass-shell while the corresponding partons 
in the parton shower might be virtual particles who branches into new particles. 
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Figure 9: distribution of ^gnai — ^^forced • 


The difference is fortunately small, roughly 10“®, compared with the upper limit 
Zmax — 1- Thus the kinematics is slightly different in the two cases. The difference 
in kinematics is illustrated by hgure 0 that shows the distribution of the angle 
between the parton given by the matrix elements and the corresponding forced 
parton. The angular differences are somewhat larger than the z differences but they 
are still small. The larger angular differences is mainly explained by the fact that 
the angle between the daughters in a branching decreases when the daughters obtain 
non-zero masses. 

All other branchings are treated as before except that the masses of the partons 
not included in the history are forced to have masses below a given threshold. The 
matrix elements generate no emission below the mass cut at 9 GeV. The hybrid 
mass threshold must therefore be at least 9 GeV in order to account for the emission 
missing from the matrix-element part. Since no 5-parton emission is generated by 
the second order matrix elements, one could allow a threshold higher than 9 GeV in 
order to account for this potential emission. However, if any such mass is larger than 
one of the forced masses, the result would be a different history than the chosen one. 
Thus one possible mass threshold is to choose it as the smallest of the two forced 
masses. Another choice is a hxed mass threshold at 13 GeV, a value that gives the 
same average multiplicity as the original parton shower. This may be viewed as a 
pragmatical compromise between the two extremes above. 

To show the exact execution the forced parton-shower algorithm is traced on the 
diagram shown in hgure H. 

a The initial partons 1 and 2 are stored in the event array. Their four-momenta 
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Figure 10: Distribution of the angle between the forced parton and the correspond¬ 
ing matrix-elements parton. 

are given in the center of momentum system and their spatial momentum is 
rotated so it coincides with the directions of the initial partons given by the 
matrix element and the chosen history. The directions are in this case given 
by Ps + P 5 + Pq and p 2 - The partons are on the mass-shell. 

b The mass of parton 1 and the 2 ; value of 1 34 are forced. The mass of parton 

2 is chosen by solving equation (0), with the extra condition that it is below 
the previously mentioned threshold. The type of branching for parton 2 is 
also decided. Because the parton is a antiquark the only available branching 
is q —> qg. If the branching parton is a gluon then one has to choose between 
branch g—i-gg or branch g—^-qq. Once the branch type is chosen the z value is 
chosen according to one of the probabilities in ( 0 ). 

c The mass and z value of parton 3 is chosen in the same way as for parton 2 while 
the mass and z value of parton 4 is forced. Then the ip angle of parton 3 and 
4 is forced around the axis of parton 1. 

d Later on, after the masses of 5 and 6 have been selected, the ip angle of parton 5 
and 6 is forced. Then the parton shower is allowed to evolve in its usual way 
until all the partons are left on their mass-shells. 


4 Results 

In order to evaluate the hybrid between matrix elements and parton showers complete 
events are generated with parton showers, matrix elements and the hybrid. All events 
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Figure 11: Distributions of charged particle multiplicity 


are generated at the Z°-resonance 91.2 GeV. As jet-finding algorithm the JADE scheme 


(technically the PO variant |M) is used. Only four-jet events are kept for further analysis. 


The three different methods have been compared in the following ways: 


1. The multiplicity distributions of charged particles are shown in figure 0- The hy¬ 
brid is executed with three different mass thresholds: fixed at 13 GeV (ME-I-PS I), 
fixed at 9 GeV (ME-I-PS II) and minimum of the two forced masses (ME-I-PS III). 
The lower limit on jet separation is set to djoin = 0.02 . This choice of djoin give well 
separated jets and reasonable high event statistics. The distribution of the parton 
shower (PS) and the three hybrid distributions nearly coincide. These four distribu¬ 
tions are also in good agreement with the experimental data [|T^] from the DELPHI 
detector at the LEP collider. The only distribution that is apparently different and 
doesn’t agree with the experimental data is the matrix elements distribution (ME). 


2. The subjet multiplicity is plotted in figure 0 - The number of jets is plotted as a 
function of the lower limit on the jet separation djom < 0.03 for events with four 
jets at djoin = 0.03 . Thus the subjet multiplicity gives a picture of how narrow 
the particles are in the four original jets. The events have been generated using the 
fixed mass threshold at 13 GeV. The curve corresponding to the matrix elements 
has the strongest tendency towards 4 jets. This is due to the fact that the matrix 
elements only generates 4 partons; further jets can only appear by fluctuations in 
the fragmentation stage. 


3. The angular distributions defined in equations (0-0) are plotted in figures 0 to 
0. The events have been generated using the fixed mass threshold at 13 GeV and 
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Figure 12: subjet multiplicities 


with the lower limit on jet separation set to djoin = 0.03 . The distributions of the 
matrix element (ME) and the hybrid (ME+PS) nearly coincides for all the three 
angular variables while the distribution of the parton shower (PS) clearly deviates. 
Comparison with data is complicated, since it is difficult exactly to reproduce the 
experimental selection criteria, but we know fl^ that the matrix elements describe 
them well. Thus |cos0Jjj^| distributions given by the hybrid, the parton shower 


and by experimental data from the OPAL detector at LEP are scaled by the 
distribution given by the matrix elements. The experimental data is scaled by the 
ME distribution given in the referred paper while the other two distributions are 
scaled by the ME distribution shown in figure |l^. The result is shown in figure 
The hybrid distribution clearly matches the experimental distribution better than 
the parton-shower distribution does. 


5 Summary and Outlook 

The results above shows that the hybrid simulation of hadronic e+e“ events gives almost 
the same results as the matrix elements in the case of the three angular variables cos+bz, 
Icos^^pj^l and cosa; 34 . The hybrid also matches the parton showers closely when it comes 
to the charged particle multiplicity and the subjet multiplicity. This indicates that the 
hybrid gives the best total description of the hadronic e’''e“ four-jet events at the Z° 
energy resonance. 

A possible continuation of the work described in this paper is to expand the hybrid 
to simulations of the full second order event qq+qqg+qqgg+qqqq. In such events the 


16 













ME+PS 

ME 

PS 


2.4 

2.2 

2 

1.8 

1.6 

1.4 

1.2 

1 

0.8 

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

COS XBZ 

Figure 13: cosxbz distributions 
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Figure 14: |cos6*nr| distributions 
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Figure 15: 0030:34 distributions 
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Figure 16: |cos 6 *J^j^| distributions scaled by the ME distribution 










scale ambiguity in Og affects the rate of 4-parton events while the shape of these 4-parton 
events are unaffected. The hrst order matrix elements correction to parton showers has 
already been investigated by Michael H. Seymour |^. Another much more difficult task 
is to apply the forced parton showers method to the description of hadron collisions. 
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